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Preface 


The  multi-faceted  mission  of  the  Space  Shuttle  necessitates 
the  orbiter's  ability  to  change  from  one  orbit  to  another. 

Many  times  an  orbit  change  includes  a  change  in  orbit  inclin¬ 
ation  which  is  usually  accomplished  by  a  rocket  burn  in  space. 
The  Space  Shuttle  Orbiter  can  possibly  achieve  changes  in 
orbit  inclination  by  making  use  of  the  orbiter's  airplane¬ 
like  ability  in  the  upper  atmosphere.  This  study  examines 
this  possibility  and  compares  it  with  the  rocket  maneuver  in 
space. 

Most  of  the  time  spent  on  this  study  involved  the  opti¬ 
mization  of  such  an  aerodynamic  maneuver.  Unfortunately,  no 
results  from  the  optimizat ion  problem  ensued;  however,  the 
problem  is  discussed  and  a  better  problem  statement  is  proposed 
in  Chapter  VII. 

Thanks  are  due  to  Captains  William  Wiesel  and  James  Rader 
for  their  advice  concerning  the  system  dynamics  and  the  optimum 
control  problem.  Captain  Ray  Barker  was  helpful  in  establishing 
a  model  atmosphere.  Special  thanks  go  to  my  wife,  Wendy,  for 
her  expertise  in  editing  and  typing  the  final  report. 


Roger  R.  J.  Harding 
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Abstract 

Two  types  of  skip  reentry  trajectories  are  examined  in 
the  70  km  to  95  km  altitude  region.  The  first  is  a  maximum 
lift-to-drag  analysis  which  indicates  that  an  aerodynamic 
maneuver  in  order  to  change  the  orbit  inclination  is  profitable 
when  compared  to  a  rocket  burn  in  space  to  effect  the  same 
change  in  orbit  inclination.  The  maximum  changes  in  orbit 
inclination  achieved  aerodynamical ly  were  approximately 
.8  degrees.  The  second  type  of  analysis  considered  the  optimal 
control  problem  for  a  skip  reentry  trajectory.  The  specific 
problem  posed  was:  Find  the  optimum  angle  of  attack  and 
bank  angle  controls  which  minimize  the  amount  of  work  done 
by  drag  for  a  specific  change  in  orbit  inclination.  No  results 
were  obtained  from  this  analysis  due  to  the  problems  encountered 
when  the  optimization  technique  was  applied  to  the  specific 
problem. 
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OPTIMUM  ORB II  PLANE  CHANGE  USINu  A  SKIP  REENTRY  TRAJECTORY 
FOR  THE  SI’ ACE  SHUTTLE  OR  BITER 


1  Introduction 


Foreword 

The  Space  Shuttle  is  a  spacecraft  system  comprising 
three  main  elements;  two  solid  propellant  rocket  boosters, 
an  external  tank,  and  the  orbiter.  The  solid  propellant 
rocket  boosters  are  expended  during  the  boost  phase  of  the 
mission  and  are  recovered  foam  the  ocean  for  reuse.  The 
external  tank  remains  coupled  with  the  orbiter  until  after 
orbit  insertion;  when  the  fuel  in  the  external  tank  is 
exhausted,  the  tank  is  jettisoned  and  reenters  the  earth’s 
atmosphere.  The  orbiter  is  capable  of  delivering  a  variety 
of  payloads  into  orbit  as  well  as  retrieving  payloads  for 
return  to  earth.  The  orbiter  can  also  carry  supplies  and 
provide  accomodations  for  up  to  four  payload  specialists 
as  well  as  the  crew  of  three  (command  pilot,  copilot  and 
mission  specialist).  Upon  completion  of  orbital  operations 
the  orbiter  reenters  earth's  atmosphere  and  terminates  the 
mission  with  an  airplane- 1 ike  landing. 

The  purpose  of  this  study  is  to  examine  the  possibility 
of  using  the  orbiter' s  aerodynamic  capabilities  to  change 
the  orbit  inclination.  Previously,  spacecraft  have  changed 
orbit  inclination  by  thrusting  with  rockets  perpendicular  to 
the  orbital  plane.  Such  a  burn  is  expensive  in  terms  of 


fuel  needed  for  the  mission  and  therefore  also  expensive 
in  terms  of  added  payload  weight  during  the  boost  phase. 

The  orbiter  is  unique  in  that  it  has  a  relatively  good 
hypersonic  lift-to-drag  ratio  (L/D)  of  approximately  2  (see 
Ref  8  ).  A  high  L/D  indicates  a  possibility  that  an  aero¬ 
dynamic  maneuver  during  a  skip  reentry  trajectory  could 
provide  an  effective  means  for  changing  the  orbit  inclination. 

A  skip  reentry  maneuver  involves  entering  the  earth’s 
atmosphere,  performing  the  desired  aerodynamic  maneuver  and 
"skipping,"  back  out  of  the  atmosphere. 

back-of-the-Envelope  Calculat ion 

The  relatively  high  L./D  for  the  orbiter  is  the  motivation 
for  this  study.  A  simplified  calculation  of  an  orbit  plane 
change  using  L/D  equal  to  2  will  follow;  this  calculation 
assumes  all  motion  lies  in  a  plane  and  the  changes  in 
velocity  due  to  the  aerodynamic  forces  are  instantaneous. 

A  typical  reentry  velocity  is  7200  meters/second;  a  loss  of 
300  meters/second  is  assumed  to  be  a  maximum  acceptable,  since 
300  meters/second  is  the  velocity  deficit  which  can  be  negated 
using  the  orbiter's  maneuvering  system  (see  Ref  8  ) . 

Example 

vinitial  =  *  *  7200  ™/s  ’•  vfinal  =  7200  '  300  =  6900  m/s 

y 


initial 


0 


X 


Using  Newton’s  Second  Law,  acceleration  =  (  Forces) /mass 

*x’  =  Ax  =  300  ru/s  =  drag/orbiter  mass 
y  =  Ay  =  lift/orbiter  mass 

Using  these  two  equations, 

lift/drag  =  2  =  Ay/ (  300  m/s) 
or 

Ay-  600  m/s 

From  the  geometry  of  the  problem, 

total  A  V  =  (6002  +  3002 =  670  m/s 
so 

i  =  sin* *  A  y/Vp 
1  25  5° 

The  velocity  of  the  orbiter  lies  in  the  orbit  plane, 
thus  the  change  in  direction  of  the  velocity  vector,  fyi,  is 
also  the  change  in  the  orbit  plane.  A  plane  change  of 
approximately  5  degrees  can  be  expected  from  such  a  skip 
turn  maneuver  according  to  this  very  rough  analysis. 
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Summary 

This  study  examines  a  skip  turn  reentry  maneuver  in  order 
to  effect  a  change  in  orbit  inclination  for  the  Space  Shuttle 
Orbiter. 

Second-order  parameter  optimization  finds  optimum  time 

histories  for  the  controls,  which  in  this  case  are  angle  of 

attack  and  bank  angle.  The  optimization  scheme  finds  these 

controls  as  polynomial  functions  of  time  while  minimizing  a 

performance  index  (work  done  by  air  drag) .  The  problem  is 

further  constrained  by  specified  final  conditions,  which  in 

this  case  are  maximum  orbital  plane  change  and  attainment 

of  a  specific  altitude  upon  completion  of  the  maneuver.  The 

theory  for  second-order  parameter  optimization,  or  suboptimal 

control  as  the  procedure  is  often  called,  is  presented  in 
Chapter  VI. 


Before  suboptimal  trajectories  can  be  determined,  it  is 
necessary  to  have  a  model  atmosphere,  the  orbiter  equations 
of  motion,  and  the  orbiter' s  aerodynamic  data. 

The  conditions  of  the  earth’s  upper  atmosphere  must  be 
suitably  modelled  in  order  to  evaluate  the  aerodynamic  forces 
of  lift  and  drag  acting  on  the  orbiter.  Lift  and  drag  are 
both  functions  of  atmospheric  density,  temperature,  kinetic 
viscosity  and  molecular  weight;  the  1962  Standard  Atmosphere 
provides  reasonable  approximations  of  same.  The  orbiter 
aerodynamic  data  obtained  from  NASA  is  based  on  the  '62 
atmosphere;  hence  the  1962  Standard  Atmosphere  is  the  model 


Solution  of  the  optimal  control  problem  requires  the 
integration  of  the  system  equations  of  motion  from  a  prescribed 
initial  state  to  a  final  state.  The  equations  of  motion  are 
derived  by  transforming  the  forces  acting  on  the  orbiter 
from  a  body-fixed  reference  frame  to  the  geocentric,  inertial 
reference  frame  and  then  applying  Newton's  Second  Law.  The 
forces  acting  on  the  orbiter  are  gravity,  lift,  and  drag. 

The  equations  of  motion  are  thus  dependent  on  the  model 
atmosphere. 

The  orbiter  aerodynamic  data  necessary  for  this  study 
are  the  lift  and  drag  coefficients.  The  coefficients  are 
dependent  upon  the  atmospheric  properties,  the  orbiter* s 
aerodynamic  configuration  (angle  of  attack),  and  the  state  of 
the  earth-orbiter  system. 

These  three  prerequisites  necessary  for  the  optimization 
scheme  are  quite  different  when  taken  individually  but  become 
interdependent  when  applied  to  this  problem.  The  equations 
of  motion  are  dependent  on  the  aerodynamic  forces  which  in 
turn  are  dependent  on  the  orbiter' s  aerodynamic  coefficients 
and  the  model  atmosphere.  The  orbiter  aerodynamic  coefficients 
are  also  dependent  on  the  model  atmosphere  and  the  state  of 
the  system.  The  state  of  the  system  is  evaluated  by  integrating 
the  equations  of  motion.  The  integration  is  begun  by  choosing 
an  initial  state  and  allowing  a  numerical  integration  scheme 
to  propagate  this  initial  state  forward  in  time. 

Suboptimal  control  requires  integrating  the  equations 
of  motion  using  guessed  polynomials  for  the  controls.  The 


performance  index  and  final  state  are  recorded.  The  coefficients 
in  the  control  polynomials  are  changed  slightly  and  the  equations 
of  motion  are  integrated  again  using  the  new  control  polynomials. 
Central  differencing  uses  changes  in  both  performance  index  and 
final  end  conditions  with  respect  to  the  change  in  polynomial 
coefficients  in  order  to  evaluate  the  first  and  second  partial 
derivatives.  These  partial  derivatives  are  used  in  suboptimal 
control  theory  to  make  a  better  guess  at  the  coefficients  in 
the  control  polynomials.  The  process  is  repeated  until  changes 
in  the  control  polynomial  coefficients  become  very  small. 

For  a  prescribed  final  altitude  and  orbit  plane  change* 
optimum  control  polynomials  are  found  which  yield  the  minimum 
work  done  by  drag  on  the  orbiter. 
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II  Model  Atmosphere 


Background 

The  aerodynamic  forces  action  on  the  Space  Shuttle 
Orbiter  are  a  function  of  the  state  of  the  system,  i.e. 
altitude  and  velocity,  knowledge  of  atmospheric  properties 
at  low  density  altitudes  is  essential  for  hypersonic  aero¬ 
dynamic  trajectory  calculations.  The  1962  Standard  Atmosphere 
used  by  this  study  provides  the  necessary  atmosphere  properties 
The  aerodynamic  forces  are  evaluated  by  solving  the  following 
equations . 

Lift  =  i2  V/»CLS 
Drag  =  %  V/*C0S 

where 

P  =  local  atmospheric  density 

V  =  free  stream  velocity  relative  to  the  orbiter 
Ct=  lift  coefficient 
C0=  drag  coefficient 
S  =  reference  area 

Lift  and  drag  are  dependent  on  atmospheric  properties 
in  two  ways:  (1)  directly  dependent  on  air  density,  and 
(2)  indirectly  dependent  on  the  temperature,  molecular 
weight,  density  and  kinetic  viscosity  used  for  evaluating 
lift  and  drag  coefficients.  Thus  an  accurate  model  atmosphere 
is  essential  in  order  to  evaluate  the  lift  and  drag  forces 
and  the  resulting  trajectory. 


The  1962  Standard  Atmosphere 

The  1962  Standard  Atmosphere  model  is  an  updated  version 
of  earlier  atmospheres  taking  advantage  of  increased  knowledge 
and  more  accurate  determinations  of  basic  quantities.  This 
model  extends  to  700  km  altitude  and  assumes  an  ideal  gas 
devoid  of  moisture,  water  vapor,  dust  and  obeys  the  perfect 
gas  law.  Further,  the  *62  model  is  idealized  to  a  mid-latitude 
year-round  mean  over  the  range  of  maximum  and  minimum  solar 
activity  (Ref  9  ).  Temperature  variations  and  molecular  weight 
variations  with  altitude  were  obtained  experimentally. 

The  1962  Standard  Atmosphere  is  used  in  this  study  for 
two  reasons!  the  Space  Shuttle  Orbiter  aerodynamic  data  is 
based  on  the  *62  atmosphere,  and  this  model  atmosphere  is 
accurate  and  easy  to  implement  in  this  study. 

Defining  Equations 

Atmospheric  parameters  of  interest  in  this  study  are 
temperature,  molecular  weight  of  air,  density  and  viscosity. 
Upper  atmosphere  winds  are  also  of  interest;  however,  due  to 
the  complexity  of  the  winds  (variations  with  time  and  position) 
the  model  simplifies  the  winds  to  an  atmosphere  rotating  with 
the  earth. 

The  variation  of  temperature,  T,  with  altitude  is  taken 
to  be  a  piecewise  continuous  function  with  constant  gradients, 
Lm.  This  model  atmosphere  uses  molecular  temperature,  Tm,  as 
the  defining  property; 

Tm  *  T(”°) 

Molecular  weight  of  air,  M,  is  a  constant,  M0,  up  to  90  km 
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altitude.  For  this  altitude  regime,  the  molecular  temperature 
and  the  temperature  are  equal.  The  interpolating  function 
for  the  molecular  temperature  isi 


Tm  ~  Tmb  +  Lm  -  ^b) 

The  subscript  b  denotes  a  base  altitude  which  is  the  altitude 
where  the  local  temperature  gradient  changes. 

Values  for  the  molecular  weight  of  air  are  available 
from  the  model  atmosphere  only  at  the  base  altitudes.  However, 
the  molecular  temperature  accounts  for  the  varying  molecular 
weight.  All  parameters  involving  molecular  weight  can  be 
expressed  in  terms  of  temperature  and  molecular  temperature. 

The  derivation  of  the  density  equation  uses  the  perfect 
gas  law,  the  standard  barometric  equation,  and  geopotential 
altitude.  Due  to  the  variation  of  gravitational  acceleration 
with  altitude,  geopotential  altitude  represents  a  modified 
altitude  to  account  for  this  variation.  This  study  assumes  that 
the  geopotential  altitude  is  equal  to  the  geometric  altitude. 
This  simplifying  assumption  results  in  a  negligible  five  to 
seven  percent  difference  from  the  1962  Standard  Atmosphere 
at  high  altitudes. 

Two  equations  for  density,/^,  are  derived i 

yr  exp  -  Q  i?  ~  Zbl  for  g  =  o 

ro  ro  Tmb  m 

and 

“2T  =  *3^  exP  0  +  £-)  (log-1  (^r-b))  for  Lm  \  0 


where 


/»o  =  atmospheric  density  at  sea  level 
//>b//o  =  known  ratio  at  the  base  altitudes 
Q  =  (>V1)  /  ( r2R) 

The  scale  height,  Q,  is  dependent  on  the  local  gravitational 
acceleration,  M%l r^,  the  molecular  weight  of  air,  M,  and  the 
universal  gas  constant,  R. 

Viscosity  is  a  necessary  parameter  for  the  calculation  of 
the  viscous  parameter,  VBAR.  Viscosity  is  assumed  to  exist 
at  all  altitudes  and  is  determined  as  a  function  of 
temperature,  Ts 

t3/2  ^ 

x  -  if  +"  s )  & 

where 

S  =  110.4  k 

/3  =  1  . 458E+06  kg./mK^s 

At  any  given  altitude  between  sea  level  and  700  km 
temperature,  molecular  weight,  density  and  viscosity  can  be 
determined  using  the  1962  Standard  Atmosphere.  These  parameters 
are  used  to  evaluate  the  aerodynamic  forces  and  to  integrate 
the  system  equations  of  motion. 
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Ill  System  Dynamics 

Evaluating  Space  Shuttle  Orbiter  trajectories  and  orbit 
plane  changes  involves  integrating  the  system's  equations  of 
motion  from  a  given  initial  state  to  a  final  state.  The  equa¬ 
tions  of  motion  are  derived  by  analyzing  system  dynamics  and 
applying,  Newton's  Second  Law. 

Assumpt ions 

The  derivation  of  the  equations  of  motion  will  involve 
the  following  assumptions: 

1 .  the  earth  is  spherical 

2.  the  earth  is  inertial 

3.  the  orbiter  is  a  point  mass  with  lift  and  drag 

acting  through  the  center  of  gravity 

4.  no  sideslip  occurs 

5.  atmospheric  wind  is  taken  to  be  the  rotational 

velocity  of  the  earth  at  orbiter's  altitude. 

These  assumptions  greatly  simplify  the  dynamics  without 
detracting  substantially  from  the  credibility  of  the  study. 
The  assumptions  that  the  earth  is  spherical  and  inertial  are 
reasonable  approximations  since  the  skip  reentry  trajectory 
maneuver  takes  a  relatively  short  timej  the  time  span  for 
such  a  maneuver  is  27  minutes  or  less.  This  time  span  is  not 
long  enough  for  the  effects  of  an  actual  non- inertial ,  non- 
sperical  earth  to  produce  an  appreciable  difference  from  the 
simplified  two-body  case  being  studied  here. 

The  assumption  that  the  orbiter  is  a  point  mass  with  the 
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aerodynamic  forces  acting  through  the  center  of  mass  is  realistic 
in  that  the  orbiter  would  be  trimmed  (i.e.  no  pitching,  yawing 
or  rolling  moments)  at  all  times  during  a  skip  reentry  maneuver. 

It  is  also  realistic  to  assume  that  no  sideslip  occurs. 

The  hypersonic  velocities  severely  limit  the  configurations 
of  an  orbital  vehicle  reentering,  the  atmosphere.  This  problem 
is  due  to  extremely  high  aerodynamic  heating,  faced  by  any  reentry 
vehicle.  Aerodynamic  heating  is  highest  on  surfaces  of  the 
vehicle  exposed  to  the  freestream  flow.  The  orbiter's  primary 
heat -protective  surfaces  are  around  the  nose  and  along  the 
underside  of  the  vehicle  with  secondary  heat-protective  surfaces 
along  the  sides  arid  top.  Thus  the  relative  velocity  vector 
must  lie  along  or  as  close  as  possible  to  the  longitudinal 
center  line  of  the  orbiter. 

Trying  to  accurately  model  the  upper  atmospheric  winds 
would  be  extremely  complex.  Such  a  wind  model  would  involve 
time  of  day,  seasons,  position,  etc.  Incorporating,  such  a 
model  goes  beyond  the  scope  of  this  study;  however,  assuming 
an  inertial  atmosphere  is  probably  not  justifiable.  Assuming 
an  atmosphere  rotating  with  the  earth  is  a  sufficient  wind 
model  for  this  study. 

Equations  of  Motion 

The  equations  of  motion  are  derived  by  establishing  the 

appropriate  coordinate  systems  and  applying  Newton's  Second 

law.  Two  coordinate  frames  are  used.  The  XYZ  frame  is  inertial 

and  fixed  at  the  earth's  center.  The  V  .ftlT  frame  is  fixed  at 

rel 

the  orbiter's  center  of  mass  (Tig  111.1). 
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The  state  of  the  system  can  be  described  by  the  state 


variables  x,y,z,and  x.y.z.  These  quantities  are  evaluated 
with  respect  to  the  inertial  frame. 

The  Vre^ML  frame  is  defined  in  terms  of  the  state  variables. 

Vrei  -  V  +  S3  X  r 


where 


L  =  M  X  V 


rel 


V  =  xi  +  yj*  +  zk  (m/s) 

_  A  A  A 

r  =  xi  +  yj  +  zk  (m) 

Z o  =  7.292115856  x  10'^  k  (radians/sec) 

A  A  A 

Unit  vectors  are  denoted  by  the  hat  symbol,  ,  and  i,j,k 
represent  unit  vectors  for  the  KYZ  frame.  The  angular  velocity 
of  the  earth  is  constant  and  denoted  by  u> .  Thus  V re ^  is  the 
velocity  of  the  orbiter  relative  to  the  atmosphere,  M  estab¬ 
lishes  a  local  horizontal  and  T  completes  the  right-hand 
coordinate  frame. 

Aerodynamic  forces  are  easily  expressed  in  the  body- 
centered  frame  (Fig  III.l). 

Lift  =  X(cos  ^  L  -  sin^  M) 

5F3J  -  -fivrel 


where 


X  =  magnitude  of  the  lift  vector¬ 
s’  ~  magnitude  of  the  drag  vector 
f  =  bank  angle  of  the  orbiter 
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The  unit  vectors  can  be  expressed  in  the  XYZ  coordinate 
frame.  A  transformation  matrix  is  needed  to  express  the  forces 
acting  on  the  orbiter  in  terms  of  the  inertial  frame.  The 

/VA 

unit  vectors,  V  jMl ,  provide  a  means  for  transforming  an 
aerodynamic  vector  in  the  body  frame  to  the  inertial  frame 
because  the  unit  vectors  are  expressed  in  the  inertial  frame. 
Thus  a  transformation  matrix,  T,  can  be  formed  by  using  the 
unit  vectors  as  column  vectors. 


T  = 


i  ■ 

A  |  A  I  A 

V  ,  I  M  '  L 
rel,  i 

i  i 


Expanding  the  defining  equations  for  the  V  ^ML  coordinate 
frame  starts  the  derivation  for  the  transformation  matrix. 

Vrel  =  (x  -  coy)i  +  (y  +  cjx)j  +  zk 

A  A  A 

=  avi  +  bv j  +  cvk 

M  =  "VrelX  r  =  ( bvz  -  cvy)  i  -  (avz  -  cvx)j  +  (avy  -  bvx)k 
=  ami  +  braj  +  cmk 

—  _  _  /v  A. 

L  =  MXVrel  =  (cvbra  -  bvcm)  i  -  (cva„,  -  avcm)j 

A 

+  (bvam  -  avbm)k 

A  A.  A 

=  axi  +  bjj  +  cxk 


Thus  the  transformation  matrix  becomes 


av/ Sy 

am/sm 

al/sl 

T  = 

> 

co 

> 

£3 

t>m/ sm 

bl/sl 

> 

CO 

> 

o 

cn/sm 

cl/si 

Si  = 

(a,2  + 

b^  +  c 

t2>% 

i  =  v, m, 1 
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where 


The  equations  of  motion  can  be  reckoned  using  Newton's 


Second  Law. 

L  F  =  nia 

where 


f  ^gravity^XYZ  +  (LltL')XYZ  +  (Dratt)xyz 


^ gravity ^XYZ  +  f(Ltft^VML  +  KDrag) 


VML 


Further , 


^gravity^XYZ 


GMm  ,  *  Z'  ^ 
— > r  (xi  -t  yj  +  zk) 
r 


where 


GM  =  earth's  gravitational  parameter 


m  =  mass  of  the  orbiter 


...  ..A  ..A  ..A 

a  -  xi  +  yj  +  zk 


(1  if t) 


VML 


0 


-  X sin  $ 
X  cos 


(Drag) 


VM1, 


8 

0 

0 


By  substituting,  into  Newton's  Second  Law  the  equations  of 
motion  can  be  found.  In  component  form  the  equations  ares 


3  ,  X  i  al  ,  "m 


"  ,  <*.  ,  J,  1-  /  1,1  \  25  “v 

x  =  -  (uM/r  )x  4-  -  (cos  <p  —  -  sinf>  —  )  -  —  — 

m  '  't;  ’  c  '  me;.. 


m 


’m 


m  Sy 


•  •  1  jZ  D-.  ^m 

y  '  -  (GM/r  )y  +  -  (cosf  —  -  sin <p  —  ) 

m  Sj  sm 


rn  Sy 


••  t  ya  t  3\  .  X  .  /  cl  ,  cm  cv 

z  =  -  ('jM/r  )z  +  —  ( cos <t>  —  -  sin <j>  —  )  -  — 

m  r  s,  r  s_  '  m  sv 


’m 
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These  three  second-order  equations  can  be  broken  down  to  six 
first-order  equations  by  redefining  the  state  variables. 


w,  = 


w-,  = 


w-,  = 


w,  = 


w  c  = 


W.  = 


The  equations  of  motion  become 


W1 

w-, 


=  W, 


=  W  r 


Wo  =  w . 


w,  = 


(UM/r^)W|  +  ^  (cos 


W'c 


-  sin 


-  (OM/r3)w2  +  *  (cos 


=  -  (GM/r3)w^  +  —  (cos 


am 

sm 


-  sm  —  ) 


m 

sm 


.5!! 

m  Sv 


g 

m  s  w 


cm  cv 

sin  —  )  -  —  — 
sra  m  sv 


The  a's,  b's  c’s  and  s's  are  also  functions  of  the  new  variables, 
These  six  first-order  equations  are  integrated  through  time 
by  the  CC6600  differential  equation  integrator,  ODE. 


Orbit  Plane  Change  Equation 

The  change  in  orbit  inclination,  or  orbit  plane  change, 
due  to  the  skip  reentry  maneuver  is  determined  by  comparing 
the  states  of  the  system  before  and  after  the  maneuver.  The 
angular  momentum  vector,  H,  involves  the  complete  state  of  the 


orbiter  and  it  is  also  normal  to  the  orbit  plane.  Thus  the 
orbit  plane  change  can  be  found  by  examining  how  the  angular 
momentum  vector  changes  direction  (Fit'  Iii.2).  The  inner 
product  for  two  vectors  is  well-suited  for  this  problem. 

Let  the  subscripts  i  and  f  denote  initial  and  final  states 
respectively. 

Ht  •  Hf  =  |TTi|  |!if|  cos  Ai 

and  indicate  the  magnitudes  of  the  initial 

and  final  angular  momentum.  Further, 

H  =  r  X  V 

A  i  =  orbit  plane  change 
Rearranging: 

.1  (riXVt)  •  (rfXVf) 

A  1  "  C°S  [riX  ^ |  |rf-xVf| 

The  plane  change  calculation  is  easily  carried  out  using  the 
above  equation. 

Integrating  the  equations  of  motion  to  a  desired  time 
and  evaluating  the  change  in  the  state  of  the  system  at  this 
time  represent  the  fundamental  dynamics  of  theproblem. 
Evaluations  of  the  aerodynamic  forces  and  the  performance 
index  are  discussed  in  the  following  chapters. 
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IV  Space  Shuttle  C rbi ter  Aerodynamics 


The  orbiter  is  assumed  to  be  a  point  mass  with  the 
aerodynamic  forces,  lift  and  drag,  acting  through  the  center 
of  mass.  This  assumption  simplifies  the  orbiter  aerodynamics 
to  evaluating  only  lift  and  drag. 

Lift  and  drag  forces  are  dependent  upon  the  state  of  the 
system;  thus  the  forces  must  be  constantly  evaluated  during 
the  integration  of  the  equations  of  motion. 

X  =  /°(x,y.z)  CL(x.y,z.x,y,z)  Vrel2(x,y,ztx,y,z) 

B  =  i?  />  (x,y,z)  CD(x,y,z,x,y,z)  Vrel  (x,y, z,x,y, z) 

The  atmospheric  density,  /\  is  obtained  from  the  model  atmos¬ 
phere  and  Vrei  is  obtained  continuously  during  the  integration 
of  the  equations  of  motion.  This  chapter  explains  how  the 
lift  and  drag  coefficients,  CL  and  ,  are  determined. 

Orbiter  Aerodynamic  Data 

The  orbiter  aerodynamic  data  was  obtained  from  the  Space 
Shuttle  reentry  office  of  NASA  in  Houston,  Texas.  The  infor¬ 
mation  included  lift  and  drag  coefficient  tables  as  functions 
of  angle  of  attack,  <* ,  and  the  viscous  parameter,  VBAR  (see 
Table  IV.l).  Also  included  in  the  data  package  was  the  infor¬ 
mation  necessary  to  evaluate  VBAR  and  an  explanation  of  how 
the  lift  and  drag  coefficients  for  the  orbiter  were  derived. 

The  viscous  parameter  is  a  subscale-fullscale  simulation 
parameter  based  on  the  freestream  properties  Mach  number,  M^, 
Reynolds  number,  Re,  and  Temperature,  T.  VBAR  is  defined  as 


(Ref  4) 


VBAR  =  MJC^Re)^ 


where 


c. 


(r’/Dk 


T  +  122.1  x  10~^5/I) 

I  +  122.1  x  10‘(5/T,) 


with 

T'  =  T  (.468  +  . 532(IW/T)  +  .l95(M)/2  M  2  ) 
K  =  .5 


Tw  -  2000  F 

r  =  1.15 


C  is  the  proportionality  factor  for  the  linear  viscosity- 
temperature  relationship,  k  is  an  empirical  constant.  The 
specific  heat  ratio,  t,  is  assumed  to  be  constant  for  the 
flight  conditions  analyzed.  A  constant  wall  temperature,  Tw, 
is  also  assumed  for  the  orbiter. 

The  freestream  Mach  number  is  the  ratio  of  Vre^  tf>  the 
frees tream  speed  of  sound-,  a^,.  The  speed  of  sound  is  assumed 
to  exist  at  all  altitudes  and  can  be  evaluated  using  the 
following  equation i 

ao*  =  (  ^  1m/M0  )  2 

where 

R  =  8.31432  x  10 3  N-m/kg  k 

The  freestream  Reynolds  number  is  based  on  the  orbiter 
length,  s  =  32.77  m. 

Re  - 

kinetic  viscosity,,*,  is  assumed  to  exist  at  all  altitudes 
and  is  obtained  from  the  model  atmosphere. 


For  values  of  VBAR  greater  than  .01  and  less  than  .08 
the  lift  and  drag  coefficients  have  been  determined  experi¬ 


| 

| 
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mentally.  For  values  of  VBAR  greater  than  .08  to  its  maximum 
value  of  5.2,  the  coefficients  are  determined  analytically 
using  the  Lockheed  Empirical  Bridging  Formula  (Appendix  A). 
This  formula  bridges  the  transitional  flow  regime  between 
continuum  and  f ree-inolecular  flow  regimes. 

These  tables  of  aerodynamic  coefficients  are  based  on 
an  assumed  nominal  reentry  trajectory  and  therefore  nominal 
velocities.  Thus,  VBAR  is  a  parameter  dependent  primarily 
on  temperature  (i.e.  altitude).  Typically,  the  extreme  values 
of  VBAR,  .01  and  5.2,  correspond  to  altitudes  of  67  km  and 
190  km  respectively. 

Given  the  state  of  the  system,  VBAR  can  be  evaluated. 

With  VBAR  and  the  angle  of  attack  available,  the  corresponding, 
aerodynamic  coefficient  can  be  found.  The  data  is  tabulated 
at  discrete  values  of  angle  of  attack  and  VBAR.  A  four  point 
bivariate  interpolation  scheme  is  employed  to  evaluate  the 
coefficients  for  angles  of  attack  and  VBAR's  between  the 
discrete  data  points.  The  bivariate  interpolation  scheme 
is  a  second-order  approximation  (Appendix  3). 
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VBAR  =  .010  .020  .040  .060  .080  . 310 


VBAR  =  .920  1.650  2.500  5.050  5.200 


Analyzing  maximum  lif t-to-dray  L/D  trajectories  serves 
two  purposes:  first,  such  an  analysis  will  indicate  whether 


or  not  a  skip  turn  reentry  trajectory  is  profitable  when 
compared  to  a  rocket  maneuver  outside  the  atmosphere;  and, 
second,  the  analysis  will  also  indicate  how  larye  an  orbit 
plane  change  can  be  expected. 

It  should  be  noted  here  that  this  analysis  is  different 
from  the  optimal  control  method  explained  in  the  next  chapter 
in  that  energy  lost  to  air  dray,  is  not  minimized,  nor  is  the 
chanye  in  orbit  inclination  maximized.  Maximum  L/D  analyses 
are  indicative  of  the  performance  of  any  aerodynamic  vehicle. 

A  maximum  L/D  trajectory  implies  that  at  each  point  alony 
the  aerodynamic  maneuver,  the  orbiter  is  oriented  to  the  anyle 
of  attack  corresponding  to  the  maximum  L/D  ratio. 

Velocity  Performance  Index 

In  order  to  effect  an  orbit  plane  chariye  in  space,  a 
rocket  burn  directed  outside  of  the  orbit  plane  is  necessary. 
Thus,  yiven  an  initial  orbit  and  a  desired  plane  chanye, 
a  corresponding  change  in  velocity,  A  V,  can  be  calculated. 

The  equation  for  this  computation  is: 

AV  =  2Vi  sin  (Ai/2) 

where 

=  velocity  prior  to  the  maneuver 
Ai  -  orbit  plane  chanye 
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Figure  V'.  I  Velocity  VecLor  Geometry 


This  equation  assumes  circular  initial  and  final  orbits  with 
equal  initial  and  final  velocities. 

Using  a  skip  turn  reentry  maneuver  to  effect  an  orbit 
plane  change  results  in  energy  lost  due  to  air  drag.  Because 
the  maneuver  starts  and  ends  at  the  same  altitude,  the  energy 
lost  is  all  kinetic  energy,  i.e.  the  maneuver  results  in  a 
AV.  The  question  of  whether  or  not  the  aerodynamic  maneuver 
is  profitable  can  now  be  properly  posed:  Given  an  orbit  plane 
change,  is  the  AV  required  for  a  rocket  maneuver , A V^,  greater 
than  the AV  required  for  an  aerodynamic  maneuver,  AV^?  Or 
asking  the  same  question  differently:  For  a  given  orbit  plane 
change,  is  the  ratio  greater  than  1?  If  the  answer 

to  these  questions  is  yes,  then  an  aerodynamic  maneuver  is 
profitable. 

The  rat io  A V^/AV^  is  called  the  velocity  performance 
index.  The  maximum  L/D  analysis  compares  the  velocity  perfor¬ 
mance  indices  for  bank  angles  varying  from  0  to  90  degrees 
at  perigee  altitudes  of  95,  90,  85,  80,  75  and  70  km.  These 
perigee  altitudes  correspond  to  initial  conditions  neglect ing 
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the  atmosphere.  Actual  perigee  altitudes  will  vary  from 
those  listed  above  due  to  atmospheric  effects. 


Discussion  of  Results 

lwo  different  types  of  graphs  are  shown  in  the  following 
figures*  velocity  performance  index  vs.  bank  angle,  and 
orbit  inclination  change  vs.  bank  angle. 

The  velocity  performance  index  for  perig.ee  altitudes 
less  than  and  equal  to  8 j  km  exceeds  one,  as  shown  in  Figures 
V .4  through  V.7.  Ihe  bank  angle  at  which  the  velocity  perfor¬ 
mance  index  becomes  greater  than  one  represents  the  minimum 
bank  angle  for  the  corresponding,  perigee  altitude  at  which 
the  aerodynamic  maneuver  becomes  profitable  over  the  rocket 
burn  maneuver  in  space. 

For  a  perig.ee  altitude  of  7  5  km,  the  praph  ends  at  70 
degrees  bank  angle  tFi,  \.b>.  If  the  bank  angle  exceeds 
70  degrees,  the  aerudviiam ic  maneuver  is  no  longer  a  skip 
maneuver;  it  becomes  a  total  reentry  where  the  trajectory 
does  not  leave  the  atmosphere.  A  similar  situation  occurs 
for  bank  angles  greater  than  50  decrees  for  a  perigee  altitude 
equal  to  70  km.  Altitudes  loss  than  70  km  are  not  within  the 
altitude  regime  for  which  orbiter  aerodynamic  data  is  available. 

Changes  in  orbit  inclination  are  insignificant  for  perigee 
altitudes  greater  than  80  km  (Fig  V.8  through  V.13).  Maximum 
obtainable  orbit  inclination  changes  for  perigee  altitudes 
less  than  80  km  are  approximately  .8  degrees.  Velocity  lost 
due  to  air  drag  is  another  parameter  of  concern.  Typically, 
velocity  lost  to  drag,  while  obtaining  .8  degrees  orbit 
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inclination  chance  ran.es  from  60  m/s  to  80  m/s.  If  800  m/s 
is  an  acceptable  velocity  loss,  then  orbit  inclination  changes 
greater  than  1  decree  are  obtainable..  This  can  be  accomplished 
by  orienting  the  initial  orbit  such  that  the  orbiter  enters 
the  atmosphere  earlier  and  leaves  later.  An  alternative  would 
be  to  perform  the  desired  orbit  inclination  chance  over  a 
series  of  skip  turn  reentry  trajectories,  i.e.  a  multiple 
orbit  revolution  maneuver. 


bank  Angle  (decrees) 


Figure  V. 


Velocity  Performance  Index  lor  95  km  Perigee 
A1 citude 


Inclination  Change  (degrees) 


Figure  V.10.  Inclination  Change  for  85  km  Perigee  Altitude 
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VI  (  [>t  Lnial  Control  Technique 


Obtaining  a  change  in  orbit  inclination  either  an  aero¬ 
dynamic  turn  or  a  rocket  burn  in  space  results  in  a  penalty. 

In  both  types  of  maneuvers  the  penalty  is  energy  that  must 
be  supplied  to  the  system  in  order  to  obtain  the  desired  inclin¬ 
ation  change  without  changing  the  other  orbit  parameters.  The 
most  efficient  method  to  change  an  orbit  inclination  outside 
the  atmosphere  is  known i  fire  the  rockets  in  a  direction 
perpendicular  to  the  orbit  plane.  On  the  other  hand,  the 
optimum  angle  of  attack  and  bank  angle  control  histories  for 
an  aerodynamic  turn  to  effect  a  change  in  orbit  inclination 
are  not  easy  to  determine  and  are  presently  unknown. 

An  aerodynamic  skip  turn  reentry  maneuver  results  in 
energy  lost  due  to  air  drag,.  In  order  to  change  only  the 
orbit  inclination,  the  energy  lost  to  air  drag  must  be  made 
up.  This  is  accomplished  by  accelerating  the  orbiter  in  the 
new  orbit  plane.  Thus,  an  optimum  aerodynamic  maneuver  requires 
that  the  energy  lost  to  air  drag  be  minimized  for  a  specified 
change  in  orbit  inclination. 

Suboptimal  Control 

Suboptimal  control  theory,  also  called  second-order 
parameter-optimization,  is  a  second-order  method  which  assumes 
that  optimal  controls  can  be  approximated  by  polynomial  functions 
of  time  (Ref  5,  6).  For  this  study  the  control  variables  are 
angle  of  attack  and  bank  angle.  The  approximating  polynomials 


are  the  Chebyshev  polynomials  which  are  chosen  because  they 
are  orthogonal  on  the  norid imens ional  time  interval,  zero  to 
one.  ihe  Chebyshev  polynomials  will  insure  that  the  matrices 
in  the  subopt imal  control  scheme  do  not  become  singular  due 
to  possible  linear  dependence  of  the  polynomials. 

In  essence,  the  subopt  imal.  control  problem  is  to  find 
the  coefficients  of  the  control  polynomials,  subject  to  the 
differential  constraints  (equations  of  motion)  and  the  pre¬ 
scribed  boundary  conditions  (initial  state,  final  altitude 
and  orbit  plane  change)  ,  which  minimize  the  performance  index. 

List  of  Symbols 

B  =  vector  of  control  polynomial  coefficients 
G  =  scalar  performance  index 
E  =  vector  of  prescribed  final  conditions 
\  =  vector  of  Lagrange  multiplier  constants 
S  =  variational  derivative 

This  study  uses  the  subscript  notation  to  denote  partial 
derivatives,  i.e.  G^  indicates  the  partial  derivatives  of 
the  scalar  performance  index  with  respect  to  the  polynomial 
coefficients . 

Subopt imal  Control  Formulation  and  Procedure 

An  augmented  performance  index,  F,  can  be  formulated 
by  combining  the  performance  index  G  with  the  equality 
constraints  E. 

F(B,A  )  =  G(B)  +  A  E(B) 
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The  Lagrange  mul tipi iers , A  ,  are  introduced  to  insure  that 
linear  dependence  does  not  occur.  The  conditions  to  be  satis¬ 
fied  at  an  extremal  point  arei 

FBr(B,A)  =  0 
E(B)  =  0 

The  numerical  procedure  starts  with  guessed  values  for  B  and 
X.  Fn  and  E  are  linearized  about  these  guessed  values. 

S  F„r  =  FbbSB  +  Eb‘sx 
S  E  =  EyS  B 

T 

These  expressions  yield  SB  and  SX  which  result  in  and  E 

moving,  toward  zero.  Further,  since  S(  )  =  (  )new  -  (  )Qld 

and  it  is  desired  that  (  )  =  0,  then 

new 

sfbt  =  -pfb‘ 

SE  =  -QE 

where  P  and  Q  are  scalar  factors  which  control  the  convergence 
process.  P  and  Q  take  on  values  from  near  zero  to  one. 
Combining,  the  previous  equations  and  solving  for  SB  and  £X 
yields : 

SX  ~(EbFbb  Eg  )  OFEgFyg  Fg  +  QE) 

SB  =  -fbb-’(i>fbt  +  EbTSX) 


where 


BB 


G 


BB 


+  A  E 


BB 


These  partial  derivatives  are  evaluated  numerically  using 
central  differencing,  techniques. 

In  summary,  the  algorithm  for  the  second-order  suboptimal 
control  scheme  is  as  follows: 

1 .  Guess  B  and  A 

2.  Integrate  the  system  equations  of  motion  from  an 
initial  state  to  a  final  state 

3.  Compute  E.  EB>  F^,  FBB 

4 .  Select  P  and  Q 


5. 


If  F^  =  0,  Fg^  is  positive  definite,  and  SB  is  small, 


the  method  has  converged 

6.  If  convergence  criteria  have  not  been  satisfied, 

set  B  =  B  +  SB  and  A  -  A  +  SA  and  iterate  procedures 
2.  through  6, 

The  procedure  is  straightforward,  but  some  discussion  is 
necessary  concerning  initial  guesses  for  A  and  B  and  selecting 
P  and  Q. 

A  good  initial  guess  for  A  can  be  obtained  from  a  gradient 
approach.  The  equations  for  the  gradient  method  are  obtained 
by  setting  FBB  =  I  and  S A=  0.  By  manipulating  the  previous 
equations,  an  expression  for  A  can  be  found: 


A  -  (EBEBr)_1  (Q/P)E  -  EBGJ}r 


The  iteration  procedure  can  be  started  with  small  values 
for  P  and  Q  (approximately  .1).  When  it  is  apparent  that  the 
procedure  is  converging,  P  and  Q  can  be  increased.  As  soon  as 
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the  end  conditions  are  satisfied,  Q  should  be  set  equal  to 
one  in  order  to  preserve  the  final  conditions. 

The  guess  at  the  coefficients  contained  in  the  B  vector 
must  be  close  to  the  actual  optimal  solution.  The  method 
used  in  this  study  is  to  initially  consider  the  case  for 
zero  plane  change  with  a  constant  angle  of  attack.  Thus 
only  one  coefficient  and  one  control  polynomial  need  be  con¬ 
sidered.  The  optimal  constant  angle  of  attack  is  a  good  initial 
g’uess  for  a  linear  angle  of  attack  polynomial.  When  the  process 
converges  on  the  optimum  linear  angle  of  attack  polynomial, 
this  resulting  polynomial  is  used  as  an  initial  guess  for 
a  quadratic  control  polynomial,  and  so  on  until  the  next 
higher  order  polynomial  does  not  substantially  change  the 
control  history.  The  next  step  involves  specifying  a  small 
change  in  orbit  inclination  and  making  a  guess  at  a  constant 
bank  angle  history  while  using  the  previous  optimal  angle  of 
attack  polynomials.  Once  the  process  converg.es  for  this  new 
B,  the  suboptimal  control  method  is  continued  with  a  constant 
bank  angle  history  while  specifying  successively  larger  changes 
in  orbit  inclination  until  the  end  conditions  cannot  be  met. 
Using  this  B  as  an  initial  guess  for  the  largest  orbit  inclin¬ 
ation  change,  higher  order  polynomials  for  the  bank  angle 
control  can  be  found  in  a  manner  similar  to  that  described 
previously  for  the  angle  of  attack  control  polynomial. 


Vll  Conclus Lons  and  Reconimendat ions 


Maximum  L./D  Analysis 

The  analysis  and  consequent  results  in  Chapter  V,  though 
approximate,  indicate  that  a  skip  turn  reentry  maneuver  is 
profitable  for  all  perigee  altitudes  between  70  and  95  km. 

The  resulting  changes  in  orbit  inclination  are  approximately 
.8  decrees!  however,  larger  changes  can  be  achieved  if  the 
trajectory  of  the  orbiter  allows  for  more  time  spent  in  the 
atmosphere. 

Suboptimal  Control  Analysis 

Most  of  the  time  spent  on  this  study  concerned  the  optimal 
control  problem  which,  unfortunately,  yielded  virtually  no 
results.  The  problem  considered  in  this  study  was:  Find  the 
angle  of  attack  and  bank  angle  control  histories  that  minimize 
the  work  done  by  air  drag  while  meet  in;’  the  specified  end 
conditions,  orbit  plane  change  and  final  altitude.  Though 
this  problem  is  certainly  relevant  and  important,  it  is  poorly 
suited  to  suboptimal  control.  Three  major  problems  occurred 
during  the  course  of  this  study.  First,  the  final  time  had 
to  be  specified.  A  good  value  for  the  time  could  be  found 
in  the  case  of  zero  plane  change:  the  time  the  orbiter  spent 
between  entering  and  exiting  the  earth’s  atmosphere.  However, 
once  an  orbit  plane  change  was  specified,  a  good  value  for  the 
time  was  no  longer  available.  Second,  the  final  altitude 
had  to  be  specified.  This  problem  came  into  focus  when  it  was 
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discovered  that  the  suboptimal  control  scheme  was  causing,  the 
orbiter  to  lose  energy  just  so  the  final  altitude  constraint 
could  be  met.  The  unsolved  question  here  is:  How  is  a  final 
altitude  established  which  is  compatible  with  the  desired 
orbit  plane  change'/ 

Third,  at  the  hi,, her  altitudes,  95  and  100  km,  the  change 

in  the  performance  index  with  respect  to  the  control  polynomial 

_2 

coefficients  was  on  the  order  of  10  when  the  coefficients 
started  approaching  the  optimal  solution.  This  led  to  diffi¬ 
culties  in  evaluating  ihe  second-order  partial  derivatives  which 
are  necessary  for  the  second-order  suboptimal  control  scheme. 

The  problem  can  be  restated  so  that  these  difficulties 
are  alleviated.  The  restated  problem  follows:  Find  the 
angle  of  attack  and  bank  angle  time  histories  that  maximize 
the  change  in  orbit  inclination  (minimize  -Zki)  for  a  given 
acceptable  loss  in  energy.  Instead  of  integrating  the  equations 
of  motion  to  a  fixed  final  time,  integrate  to  a  final  altitude 
(equal  to  the  initial  altitude).  Thus  the  loss  in  energy  due 
to  drag  would  be  all  kinetic  energy.  This  will  eliminate  the 
first  two  difficulties  encountered  with  the  current  problem. 

In  order  to  alleviate  the  third  difficulty  it  is  recommended 
that  the  analysis  be  restricted  to  90  km  perigee  altitude  and 
lower.  The  probable  difficulty  with  this  newly-proposed 
problem  will  be  establishing  typical  losses  in  energy  that 
must  be  prescribed.  This  problem  is  not  insurmountable. 

Typical  values  for  energy  loss  can  be  found  using  an  analysis 
similar  to  that  formed  in  Chapter  V  of  this  study. 


Due  to  the  importance  and  interesting  nature  of  a  study 
such  as  this  one,  it  is  recommended  that  optimum  skip  turn 
reentry  trajectories  be  investigated  further.  The  newly- 
posed  problems  stated  previously  are  relevant  and  also  warrant 
further  and  more  accurate  investigation. 
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Appendix  As  Lockheed  Bridr.ini',  Formula 


The  Lockheed  Bridging  Formula  is  used  to  evaluate  the 
aerodynamic  force  coefficients  for  the  Space  Shuttle  Crbiter. 
The  formula  bridges  the  transitional  flow  regime  from  continuum 
flow  to  free-molecular  flow.  The  formula  is: 


trans 


=  C 


cont 


+  (C 


F.M. 


"  CcontJ 


sinn(  tt(A  +  B  logj0Kn)) 


where 


C  _  =  viscous  force  coefficient  values  at  VBAR  =  0.08 
cont 

CL.  ..  =  viscous  force  coefficient  values  at  VBAR  =  5.2 

r  .M. 

(free-molecular  flow) 


n  =  2 
A  =  3/8 
B  =  1/8 

kn  =  Knudsen  number  =  A /^ref 
Lre^  =  12.059  meters 

?  - 1 

A  =  mean  free  path  =  RT/P  (  2  N<r  ) 

3 

R  =  universal  gas  constant  =  8.314  x  10  N-m/kg  k 
T  =  temperature 
P  =  pressure 

26  “1 

N  =  Avogadro's  number  =  6.022  x  10^  kmol 
O'  -  effective  molecular  collision  diameter 
=  3.65  x  10'10 


meters 


Appendix  b ;  bivariate  Interpolation 


The  four  point  bivariate  interpolation  scheme  is  used 
to  find  the  aerodynamic  coefficients  between  the  tabulated 
values.  The  interpolating  tormula  is-. 


f(x,y)  = 


where 


f(xQ  +  x,  yQ  +  y )  ~  (1  p)(l  -  q)  f(x0,y0) 
+  p(l  -  q)  f(x1,yt))  +  q(  1  -  p)  fCxQ.yj) 

+  pq  f (X| ,y| )  +  © (h2) 


=  S  x/  Ax 
=  S  y/ Ay 


The  four  point  interpolation  scheme  can  be  show  graphically 


Lf(xo*yl)  •ftxl‘yl) 


*f(x,y)=f(xo+  x,y(j+  y) 


Lf  (xo,yo)  /(Vyo} 


where 


=  Xj  -  xQ 

=  y\  -  yo 

=  x  -  Xo 


=  y  -  y« 


For  this  study,  x  and  y  correspond  to  VBAR  and  angle  of 
attack  respectively. 
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